Astronomy & Astrophysics manuscript no. d2ov8 
April 11, 2013 



©ESO2013 



Heavy water stratification in a low-mass protostar ^ 

A. Coutens''2, C. Vastel^'^, S. Cazaux^ S. Bottinelli^'^, E. Caux^-^, C. Ceccarelli'*, K. Demyki-^, V. Taquet'^, and V. 

Wakelam^ 



cn 

o 

;h 

< 
O 



< 

O 
6 



Universite de Toulouse, UPS-OMP, TRAP, Toulouse, France 

e-mail: acoutens@nbi . dk , cvastelOirap . omp . eu 

CNRS, Institut de Recherche en Astrophysique et Planetologie, 9 Av. Colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France 

Kapteyn Astronomical Institute, P.O. Box 800, 9700AV Groningen, The Netherlands 

Institut de Planetologie et d' Astrophysique de Grenoble, UMR 5274, UJF-Grenoble 1/CNRS, 38041 Grenoble, France 

CNRS and Universite de Bordeaux, Observatoire Aquitain des Sciences de I'Univers, 2 rue de I'Observatoire, B.P 89, F-33271 

Floirac, France 



Received 20 December 2012; accepted 19 March 2013 



ABSTRACT 



Context. Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular deuterium fractionation has been found 
in the environments of low-mass star- forming regions and, in particular, the Class protostar IRAS 16293-2422. 
Aims. The key program Chemical HErschel Surveys of Star forming regions (CHESS) aims at studying the molecular complexity of 
the interstellar medium. The high sensitivity and spectral resolution of the HerschelfHWl (Heterodyne Instrument for Far-Infrared) 
instrument provide a unique opportunity to observe the fundamental li.i-Oo.o transition of ortho-D20 at 607 GHz and the higher 
energy 2i 2-lo,i transition of para-D^O at 898 GHz, both of which are inaccessible from the ground. 

Methods. The ortho~D20 transition at 607 GHz was previously detected. We present in this paper the first tentative detection for the 
para-D^O transition at 898 GHz. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed 
line profiles of D2O with the same method that was used to reproduce the HDO and H2'*0 line profiles in IRAS 16293-2422. 
Results. As for HDO, the absorption component seen on the D2O lines can only be reproduced by adding an external absorbing 
layer, possibly created by the photodesorption of the ices at the edges of the molecular cloud. The D2O column density is found to be 
about 2.5 X 10'- cm"^ in this added layer, leading to a D2O/H2O ratio of about 0.5%. At a 3cr uncertainty, upper limits of 0.03% and 
0.2% are obtained for this ratio in the hot corino and the colder envelope of IRAS 16293-2422, respectively. 

Conclusions. The deuterium fractionation derived in our study suggests that the ices present in IRAS 16293-2422 formed on warm 
dust grains (~ 15 - 20 K) in dense (~ 10** - 5 x lO"* cm"') translucent clouds. These results allow us to address the earliest phases of 
star formation and the conditions in which ices form. 

Key words, astrochemistry - ISM: individual (IRAS 16293-2422) - ISM: molecules 
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1. Introduction 

Despite the low elementa l deuterium abunda nce in the local in- 
terstellar gas (1-2 X 10"^; Linsky et al.|2006 1, numerous deuter- 



ated molecules have been detected with enhanced D/H ratios in 
the environments of low-mass star-forming regions. The dou- 
bly deuterated form of water, D2O, has been in particular de- 
tected in the solar-type protostar IRAS 16293-2422 (hereafter 
IRAS 16293), which is well known for its high deuterium frac- 
tionation: the fundamental li,o-lo,i para-D20 transition at 317 
GHz was detected by |Butner et al.| ( |2007[ ) using the James 
Clerk Maxwell Telescope (JCMT), whereas the fundamental 
li,i-Oo,o ortho-D20 transition at 607 GHz was discovered with 
the Heterodyne Instrument for Far-Infrared (HIFI) instrument 
onboard the Herschel Space Observatory by Vastel et al. (2010 1. 
The observed line profile at 317 GHz shows a component in 
emission in addition to a deep and narrow absorption. The emis- 
sion component has been attributed to heavy water in the hot 
corino of this source, where the grain ices are sublimated and 
released into the gas phase ( [Ceccarelli et al.||2000| [Bottinelh 



Based on HerschelfUlFl observations. Herschel is an ESA space 



observatory with scientific instruments provided by European-led prin- 
cipal Investigator consortia and with important participation from 



NASA. 



et al.|2004) l. The absorption component, whose linewidth is 0.5 
kms~\ is seen in both transitions at the velocity of 4.2 km s"' . It 
probably originates in the foreground gas (molecular cloud and 
cold envelope), which is studied through the fundamental HDO 
transition at 465 and 894 GHz ( |Coutensetal.|2012| l. Both transi- 
tions present deep absorption features at the same velocity as the 
D2O component. The OD molecule, a key species in the forma- 
tion of deuterated water on grains, also shows a deep absorbing 
component (Parise et al. 20121 and most certainly arises from 
the same region as water and its deuterated forms. Besides, nu- 
merous absorption lines of several nitrogen molecules (NH, ND, 
NH2, and NH3) are detected towards the IRAS 16293 source 



(Bacmann et al. 2010 Hily-Blant et al. 2010), suggesting an ori- 
gin of these species in a cold environment. ,Coutens et al. ( |2012| l 



used the spherical Monte Carlo radiative transfer code RATRAN 
(Hogerheijde & van der Tak 2000) to reproduce the line pro- 
files of the two HDO fundamental transitions as well as eleven 
other detected Unes by assuming an abundance jump at 100 K 
( Fraser et al.|2001 1. From the emission line profiles, the derived 
HDO inner abundance (T > 100 K) is about 1.7 x 10"^, whereas 
the HDO outer abundance (T <100 K) is estimated at about 8 x 
10"" . A water-iich absorbing layer is, however, required in front 
of the envelope to reproduce the HDO absorption components 
observed at 465 and 894 GHz. This layer may result from the 



A. Coutens et al.: Heavy water stratification in a low-mass protostar 



photodesorption of the ices at the edges of the molecular cloud, 
as predicted by chemical models (e.g., Hollenbach et al.|2009 1. 
A non-LTE (local thermal equilibrium) modeling has been 
performed by |Vastel et al.| ( 2010| l to estimate the ortho-to-para 



D2O ratio in the cold (< 30 K) dense (< 5 x 10^ cm"^) cloud 
surrounding the IRAS 16293 protostar, adopting the density and 
temperature profiles of the envelope ( |Crimier et al.|2010l ). |Vastel] 
let al. (2010) considered the first two levels of each D2O form 



and used the computed collisional rates for the two fundamental 
deexcitation transitions of ortho- and para-D20 with para-H2 



in the 10-30 K range (Scribano et al. 2010). Their computation 



yields an ortho-to-para ratio of 1.1 ± 0.4 at a Icr level of un- 
certainty with the corresponding column densities N(ortho) - 
(8.7 + 2.1) X 10" cm-2 and N(para) = (7.8 ± 2.6) X 10'^ crn-^. 
This ratio is lower than 2.6 at a 3cr level of uncertainty, taking 
into account the statistical error as well as the overall calibration 
budget for the HIFI band lb. The comparison between the upper 
value of the measured D2O ortho-to-para ratio and the thermal 
equilibrium value shows that they are consistent with a gas at 
a temperature higher than about 15 K and, therefore, with the 
assumed absorbing gas location. 

The present paper aims at determining the D2O abun- 
dance distribution throughout the whole protostellar envelope 
of IRAS 16293, from the hot corino to the foreground water- 
rich absorbing layer (see Figure 2 in ' Coutens et al.|20I2|). Using 
the previous results obtained in a similar way by [Coutens et aT 



( |20T2| ) for HDO and H2O, the derived D2O/H2O and HDO/H2O 
ratios are then discussed and compared with a grain-surface 
chemical model (Cazaux et al. 201 1| ). The paper is organized 
as follows. First, we describe the observations as well as their 
reduction in Section 2. Then we present the modeling and the 
results in Section 3. Finally, we discuss the results in Section 4 
and conclude in Section 5. 



2. Observations 

In the framework of the Chemical HErschel Surveys of Star 
forming regions (CHESS) key program ( Ceccarelli et al.|2010 i, 
we observed the solar-type protostar 
Hersc helfHlFl instrument ( |Pilbratt et al. 
2010 1. The targeted coordinates were Q'2000 
24° 28' 34.2 



RAS 16293 with the 
20T0l|deGraauwetal.| 



^2000 



16^32'" 22^.75, 
a position at equal distance from the 



binary components A and B. A full spectral coverage of bands 
lb and 3b was performed on 2010 March 2 and 2010 March 19 
respectively, using the HIFI spectral scan double beam switch 
(DBS) mode with optimization of the continuum. The funda- 
mental ortho-D20 li,i-Oo,() transition lies in this frequency cov- 
erage at 607.349 GHz whereas the para-D20 2i,2-li,o transi- 
tion lies at 897.947 GHz. The HIFI Wide Band Spectrometer 
(WBS) was used, providing a spectral resolution of 1.1 MHz 
(-0.55 kms' at 600 GHz and -0.37 kms' at 900 GHz) over 
an instantaneous bandwidth of 4 x 1 GHz. The ortho-D20 lij- 
Oo,o transition was re-observed in February 2011 with the High 
Resolution Spectrometer (HRS) to benefit from a higher spec- 
tral resolution (~0.06 kms"^). We consequently use these new 
observations hereafter At 610 GHz (900 GHz), the beam size is 
about 35" (24" resp.) and the measured main beam efficiency is 
0.75 (0.74 resp.), whereas the forward efficiency is 0.96 in both 
cases (Roelfsema et al. 2012[ ). The DBS reference positions were 
situated approximately 3' east and west of the source. 

The data reduction of the para-D20 2i2-li,o transition at 
898 GHz was performed similarly to the HDO and H2'^0 lines 
studied by Coutens et al.| p012). This transition was observed 
eight times (four in lower sideband and four in upper sideband) 
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Fig. 1. Energy levels for the lowest rotational transitions of 
ortho- and para-D20. 



for each polarization in the selected observing mode. To produce 
the final spectra, the observations that were processed using the 
standard HIFI pipeline up to level 2 with the ESA-supported 
package Hers chel Interactive Processing Environment (HIPE) 
5.1 ( |Ott|2010| ) were exported to the GILDAS/CLASSQ software. 
We then baseline subtracted and averaged the H and V polar- 
izations at the line frequency, weighting them by the observed 
noise for each spectra. We verified for all spectra that no emis- 
sion from other species was present in the image band. This 
transition had not been discovered when |Vastel et al.| ( |2010] ) 
published their results with an early version (3.01) of the HIPE 
software, using the sideband deconvolution method on all bands. 
The present data reduction, which uses the selected spectra re- 
aligned on one single transition and averaged after an accurate 
baseline subtraction and rms weighting, leads to a lower noise 
level than the one obtained with the standard HIPE deconvolu- 
tion routine, and allows us to present the first tentative detection 
of a high-energy level transition of deuterated water (see Fig. l2]). 
For the pointed HRS data at 607 GHz, the observations of the 
H and V polarizations were processed with the ESA-supported 
package HIPE 8.0 (Ott'2010) up to level 2 and then exported to 
the GILDAS/CLASS software for the baselines to be subtracted 
and the two polarizations averaged. HIFI operates as a double 
sideband (DSB) receiver. From the in-orbit performances of the 
instrument ( Roelfsema et al.|2012 i, a sideband ratio of unity is 
assumed for the observed D2O transitions seen in absorption 
against the continuum (bands lb and 3b). The continuum val- 
ues observed over the frequency range of the whole bands are 
closely fitted by straight lines. The continuum derived from the 
polynomial fit at the considered frequency and divided by two 
to obtain the single sideband (SSB) continuum was added to the 
spectra. 

The parameters of the lines, taking into account the ortho- 
and para-D20 forms (within the CASSIS softwaraj) separated 
from the Cologne Database for Molecular Spectroscopy (Miillerl 
et al. 2005 Brunk en et aLl[2 00T), are reported in Table [T] The 
fundamental ortho-D20 transition at 607.349 GHz is clearly de- 
tected in absorption against the strong continuum (~ 0.3 K). We 
present here the first tentative detection of the para-D20 2i 2- 
lo,i transition at 897.95 GHz, seemingly in absorption against 
the continuum (~ 0.75 K). In the same table, we also report 
the parameters of the para-D20 1 i^o-lo.i fundamental line previ- 
ously observed at 316.800 GHz with the JCMT by [Butner et"aL] 
( j2007j ). We used the line decomposition in three Gaussian com- 
ponents derived by Vastel et al. (2010) : two of them are at- 



tributed to the broad emission component and the deep absorbing 

' http://www.iram.fr/IRAMFR/GILDAS 



Developed by IRAP-UPS/CNRS: http : //cassis . irap . omp . eu 
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Table 1. Derived parameters of the absorbing components of the ortho- and para-DaO fundamental lines. 



Species Transition 



Frequency 
(GHz) 



(S-' 



2.95 X 10-' 
6.29 X 10-'' 
8.61 X 10-3 



(K) 



Telescope 



(mK km/s) 



(mK) 



FWHM 
(km/s) 



Vlsr 
(km/s) 



ir.b.sl/rc 



ortho-D20 l,j- 
para-D20 li.o- 
para-D20 2] 2- 



-Oo,o 
■lo,i 
-lo,i 



607.34945 
316.79981 
897.94711 



29.1 
32.6 
60.5 



HIFI Ib/HRS 

JCMT 
HIFI 3b/WBS 



-90: 
-120: 

-734 



:49 

13 



-199: 

-220 d 
-167 d 



30 
12 



0.43: 
0.55: 

0.41 : 



0.02 
0.15 
0.04 



4.22: 
4.15: 
4.04: 



0.01 
0.04 
0.02 



0.64 ± 0.03 
0.57 ± 0.08 
0.22 ± 0.02 



line of D2O, whereas an additional component emitted at 10.1 
kms-' is likely due to the 1q-6q+ CH3OD hne at 316.795 GHz. 
This latter component was subtracted in Figurel2] The 2().2-li.i 
ortho transition at 468.25 GHz (see Fig. [T} has not been ob- 
served, and we consequently only provide a prediction with the 
modeling. 

Observations were recently carried out with the Atacama 
Pathfinder Experiment (APEX) telescope between 270 and 
330 GHz to complete The IRAS 16293-24 22 Millimeter and 
SubmilUmeter Spectral Survey (TIMASSS; |Caux et al.|[20n i. 
While the para-DaO fundamental li,o-lo,i line at 317 GHz is 
also detected in this survey, it shows a much lower signal-to- 
noise ratio than the one observed with JCMT ( |Butner et al.| 
|2007 ). This observation is nevertheless useful due to the accu- 
rate continuum at this frequency (0.235 K in T*^. The beam ef- 
ficiency is 0.73 and the forward efficiency 0.97. 



3. Modeling 

In order to determine the D2O/HDO and D2O/H2O ratios within 
the different components of this source (hot corino, colder en- 
velope, and external absorption layer), we developed the same 
modeling as performed in Coutens et al. ( |2012 i for the HDO 
and H2'^0 species with the RATRAN radiative transfer code 
(Hogerheijde & van der Tak 2000 1. We employed the ortho- 



and para-D20 colUsional coefficients calculated with ortho- and 
para-H2 by IFaure et al. ( 2012| l in the 5-100 K range. We as- 
sumed that the ortho-to-para D2O ratio is 2:1, compatible with 
the results from Vastel et al. (2010 1. Using the same parameters 
for the modeling (see details in |Coutens et al.|2012[ ), we realized 
that the continuum modeled at 317 GHz (0.38 K) is not consis- 
tent with the continuum observed by Butner et al. ( 2001) at 0.85 
K. To determine whether this divergence comes from the model- 
ing or from the observations, we used the recent APEX observa- 
tions. The SSB continuum observed with APEX is 0.31 K (T^b). 
Our modeling predicts a 0.29 K continuum value at the same 
frequency. From the comparison between transitions from over- 
lapping frequency ranges, we can estimate the calibration to be 
20%. The value estimated from the modeling is therefore con- 
sistent with the APEX continuum value, and in the following we 
will use the 0.38 K predicted by our modeling for the continuum 
of the 317 GHz transition observed with JCMT. 

With the source structure (hot corino + colder envelope) de- 
rived by jCrimier et al.| ( [2010 ), the modeling is not able to repro- 
duce any absorption feature for the D2O fundamental transitions. 
The overall absorption can only be reproduced by the external 
absorbing layer, which was discovered by jCoutens et al. ( 2012| l 
in their efforts to reproduce the absorption components seen on 
the HDO fundamental transitions. Assuming an H2 density of 
~ 10^-10"^ cm-^, a kinetic temperature of ~ 10-30 K and a b- 
Doppler parameter of 0.3 km s"^ (turbulence), the ortho-D20 
(respectively para-D20) column density in this layer is about 
1.4 X 10'^ cm-^ (respectively 1.1 x 10'^ cm"^). Our model- 
ing predicts that the 2i 2-lo,i para-D20 transition is expected 
in absorption as well, tracing the external absorbing layer The 



tentatively detected absorbing component at 898 GHz is con- 
sequently in perfect agreement with the modeling (see Fig. l2]i. 
The kinetic temperature (~ 10-30 K) has a limited influence on 
the absorbing transitions, but hydrogen density seems to have 
more influence. With a lO*" cm"^^ density, the absorption depth 
is less pronounced for the fundamental li,o-lo,i and li,i-Oo,o 
transitions at 317 and 607 GHz, requiring a higher D2O column 
density. However, the higher energy para-D20 transition absorp- 
tion is overestimated compared to the observations, suggesting 
a lower H2 density. From the above estimations, the D2O ortho- 
to-para ratio is, in the absorbing layer, about 1.3, close to the 



value determined by Vastel et al. (20101. However, taking into 
account the observation and modeling uncertainties, we cannot 
exclude that the measured value might be consistent with the 2: 1 
value at equilibrium. The quoted column densities are shghtly 
higher (less than a factor 2) than those estimated by [Vastel et al. 
(2010 1, but within their 3cr uncertainties. If this absorbing layer 
is produced by photodesorption at Ay ~ 1-4 mag as predicted by 
[Hollenbach et al.,(2009) , the D2O abundance in this layer should 
be about 6.6 x 10~^"-2.7 x 10~^, using the relation between the 
H2 column density and the visua l extiri ction, N(H2)/Av/ = 9.4 x 
10^" cm-^ mag-' (Frerking et al. 1982| l. 

To estimate the D2O abundances in the hot corino {Xm) and 
the outer envelope (Xout) of IRAS 16293, we ran a grid of models 
that take into account the structure of |Crimier et al.| P010| l as 
well as the foregound absorbing layer previously defined. The 
best-fit modeling is then obtained for an internal abundance Xi^ 
= 7 X 10-'" and an external abundance Xoui -5 x 10-'^ from 
the x^ minimization (see Fig. pj) computed on the line profiles 
according to the following formalism: 



N «chan 



X' 



_ -^ 'y^ (^obs.o ^mod.ip 



(1) 



'■=1 >=i 



with A^ the number of lines /, n^han the number of channels j for 
each line, Tobs.ij and T^odM the intensity observed and predicted 
by the RATRAN non-LTE modeling respectively in the channel 
j of the line / and nns; the rms of the line /. Here, the reduced 
X^ is equal to 2.5. At 2cr, the inner and outer abundances ranges 
are Xi„ = [9 x 10-" - 1.1 x 10-'^] and X„^, = [8 x 10-'^ - 1.1 
X 10-"]. Only upper limits can, however, be given at 3cr, with 
Xin < 1.3 X 10-'' and X^ui < 1.3 x 10-". The best-fit modehng 
is shown in Figure [2] 

To check the influence on the ortho-to-para ratio assumed 
for D2O in the protostellar envelope (i.e., the equilibrium value 
of 2: 1), we ran another grid of models with an ortho-to-para D2O 
ratio equal to 1.3 in the envelope, i.e., the same as estimated in 
the absorbing layer The results in both cases are quite in agree- 
ment and do not allow us to conclude with an ortho-to-para ratio 
different from the equilibrium value 2: 1 . In addition, we checked 
that the para-D20 emission line predicted at 317 GHz was con- 
sistent with the line observed with APEX. Regarding the 2o,2- 
li 1 ortho transition at 468 GHz, we predict that the line only 
shows an emission component with an intensity about 0.1 K 
(r^b) if observed with APEX. 
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Table 2. HDO/H2O, D2O/HDO, and D2O/H2O abundance ratios estimated towards the protostar IRAS 16293. 





Hot corino 


Outer envelope 


Photodesorption layer 




Best-fit 


3o- 


Best-fit 


30- 


Av ~ 1 - 4 mag 


HDO" 


1.8 X 10-' 


1.4- 2.4 X 10-' 


8 X 10-" 


5.5- 10.6 X 10-" 


~ 0.6 - 2.4 X 10-** 


H2O"'' 


1 X 10-^ 


4.7 - 40.0 X 10-" 


1.5 X 10-** 


7.0- 22.5 X 10-" 


~ 1.3- 5.3 X 10-' 


D,0 


7x 10-'" 


< 1.3x 10-" 


5x 10-'2 


< 1.3 X 10" 


~6.6-27x 10-1° 


HDO/H2O 


1.8% 


0.4% -5.1% 


0.5% 


0.3% -1.5% 


~ 4.8%'' 


D2O/HDO 


0.4% 


< 0.9% 


6.3% 


<23% 


~ 10.8%^ 


D2O/H2O 


0.007% 


< 0.03% 


0.03% 


< 0.2% 


~ 0.5%" 



The HDO and H2O abundances from 



Coutens et al. 



1 2012 I were recalculated with the ;jf^ minimization shown in Equation fll '' The H2O 



abundances are obtained from the H2'^0 data using the collision rates with H2 calculated by Daniel et al. J201lf. The H2'*0/H2"*0 ratio 



assumed to be 500. " The ratios remain valid if the absorbing layer is not due to the photodesorption 



0-D2O 1 i.i-Oo.o 607 GHz 
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Fig. 2. In black: D2O transitions observed with HIFI and JCMT. The CH3OD (7o-6o+) transition at 316.795074 GHz contaminating 
the para-D20 (li,o-lo,i) line profile has been subtracted for a direct comparison with the D2O best-fit model. In red: best-fit model 
obtained when adding an absorbing layer with N(para-D20) = 1.1 x lO'^ cm-^ and N(ortho-D20) = 1.4 x 10'^ cm-^. The resulting 
inner abundance is 7 x IQ-'^ and the outer abundance is 5 x IQ-'^. 
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Fig. 3. ;^'^ contours at Icr, 2cr, and 3cr obtained when adding an 
absorption layer with a para-D20 column density of 1.1 x 10'^ 
cm-^ and ortho-D20 column density of 1.4 x 10'^ cm-^. The 
best-fit model is represented by the symbol "-I-". 



4. Discussion 

4.1. Ortho-to-para D2O ratio 

At the thermal equilibrium, the ortho-to-para ratio of D2O is 
higher or equal to 2, depending on temperature (see Fig. 4 in 



Vastel et al. 2010). If the ortho-to-para D2O ratio estimated at 



1.3 in the absorbing layer really deviates from the 2:1 value as 
already suggested by | Vastel et aL] (|2010[), this could mean that 



some nuclear spin exchanges occur in the gas phase or at the 
grain surface. Indeed, the ortho-to-para ratio should be initially 
set during the molecule formation at a value reflecting the grain 
temperature (> 2). Some spin exchange mechanisms would be 
consequently required to decrease the ortho-to-para ratio. These 
mechanisms are, however, poorly known. We can imagine some 
exchange mechanisms at the grain surface with interaction be- 
tween electronic or nuclear spins. Unfortunately, the theoretical 
and laboratory measurements ( Limbach et al.|[2006 Pardanaud 
et al.||2007 1 are diflicult to carry out. In the gas phase, no fast 



radiative transition is possible between ortho and para species, 
except via proton exchange with ions. For example, to modify 
the ortho-to-para ratio of water, collisions with H^ and Ht,^ are 
required. For D2O, the following reaction could consequently 
occur via D exchange between D2O and H2D^: 



0-D2O + H2D+ 



P-D2O + H2D+, 



However, this reaction could also lead to other products for 
which the branching ratios are unknown: 



D2O + H2D+ 



HD2O+ + HD 
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D2O + H2D+ 



DiO+ 



H2 



This reaction has never been studied experimentally or theoreti- 
cally, but could be exothermic like the H2O + H3 reaction (Woon 
& Herbst|2009[l. The nuclear spin conversion could also be pos- 



sible thanks to an exchange with D"^ 



4.2. Water deuterium fractionation 



As already noticed in Coutens et al. ( 2012| l, the HDO/H2O ra- 
tios determined in the hot corino and in the external absorbing 
layer are similar (~2-5%), although the corresponding densities 
are very different, a few 10*^ cm"^ in the hot corino region and, 
typical of a molecular cloud, probably lO'^-lO'* cm'^ in the ex- 
ternal layer. Water shows a different behavior from other trac- 
ers as methanol (CH3OH) and formaldehyde (H2CO), whose 
deuteration increases with the density. Indeed, the deuteration 
of these heavier molecules that need CO ices to form is corre- 
lated with the CO depletion (hnked to the molecular hydrogen 
densit y: [Caselli et al."1999' 'Bac mann e'raL]|2002| l in prestellar 
cores ( IBacmaniTet al. 2003 , 2007||. In contrast, the water frac- 
tionation ratio does not seem to vary with density. As discussed 
in |Coutens et al. ( 2012| l, it is therefore possible that water was 
formed at low densities in the early stages of the star formation 
before the protostellar collapse. 

To understand how the deuterium fractionation of water de- 
pends on the initial conditions of the cloud in which ices are 
forming, we used in the present study a grain surface chemical 
model ( Cazaux et al.|201 1) , in which we considered the forma- 
tion of ices by accreting species from the gas phase. The icy 
mantles then grow until 99% of CO and O from the gas are de- 
pleted onto dust. We set the initial densities of O and CO in the 
gas similar to densities of O and CO in translucent clouds, i.e., 
n(0) ~ n(CO) ~ 1.5 x lO'* «« (JHollenbach et al. 2009). The 
gas phase species accrete on the dust and the chemical evolution 
of ices can be followed with time until most of O and CO are 
depleted from the gas phase (f ~ 1-5 x 10^ years, depending on 
the density of the medium). During the formation of ices, O and 
CO gas phase abundances are decreasing while the abundances 
of gas phase H2, H and D are kept constant. The density of the 
cloud sets the gas phase D/H ratio though ion-molecule chem- 
istry (deuteration of HJ sets the gas D/H ratio, as described in 
[Roberts et al.|2002] ). Therefore, a translucent cloud with an ini- 
tial density of lO-*, lO'*, and 5 x 10"^ cm"^ results in a D/H ratio 
in the gas phase of 5 X 10"^, 5 x 10"^, and 10"^ respectively 
( [Cazaux et al.|20lT ). Consequently, we assume these three ini- 
tial values for the density (10^, 10'*, and 5 x 10"* cm"^) as well as 
three different temperatures (12, 15, and 20 K). The HDO/H2O 
as well as the D2O/H2O ratios obtained during the formation of 
ices are presented in Figure p] Both ratios are strongly sensitive 
to the dust temperature. Indeed, at low temperatures, H2 is much 
more abundant than H on the dust grains, and the successive hy- 
drogenations of O with H2 dominate the formation of H2O. The 
reaction 0(^P) -1- H2 is energetically not accessible at low tem- 
peratures since the reaction has a barrier of 0.57 eV and is endo- 
ergic with an energy of 0.1 eV for H2(v=0) ( Week et al.|2006| l. 
However, this reaction becomes exo-ergic for v>0 and the barrier 
decreases to ~0.4 eV for v=l (jSultanov & Balakrishnan 2005^). 
The reaction 0(^P) + H2(v=l) proceeds through tunneling at low 
temperatures, as shown by |Weck et al. ( 2006[ l. In our model, 
we consider only the reaction of H2 vibrationally excited with 
0(^P). For a population in equilibrium at a temperature Tdus,, the 
fraction of molecules in the J - I state is ~ exp (- 150/rdust) ~ 3 
X 10"^ - 5.5 X 10""* in the range of Tdust considered in this study 



(12 - 20 K). The reaction through tunneling between H2(v=l) 
and 0(^P) dominates at very low temperatures when the dust is 
saturated with H2 molecules and the oxygen atoms present on 
the surface repetitively collide with H2. On the other hand, HDO 
forms through addition of atomic deuterium. These different for- 
mation routes imply that the HDO/H2O ratio scales with D/H2 
ratio present in the gas, while the D2O/H2O scales with (D/H2)^. 
The reaction between O and H, however, dominates at 12 K at 
short timescale because of the relatively high initial H/H2 ratio 
(~ 10"^). However, the H/H2 ratio decreases with time, since 
H atoms are quickly converted into OH and H2, and the reac- 
tion O -H H2 dominates (for f ~ 10"* years at ne = lO"* cm"-' and 
for f ~ 10^ years at «h = 10'* cm"-'). If the grain is warm, the 
H2 molecules evaporate and the formation of H2O is insured by 
successive hydrogenation with H atoms. Since the formation of 
HDO involves atomic deuterium, the HDO/H2O ratio becomes 
comparable to the D/H ratio and the D2O/H2O ratio comparable 
to (D/H)^. Therefore, the variation of the deuterium fractionation 
with the initial dust temperature is due to the different routes to 
form water molecules. In their modeling, ICazaux et al. (2011 1 
found that the deuterium fractionation of the water ice mantle 
is not altered by the cloud collapse, pointing to the fact that the 
deuterium fractionation of water must have been set during the 
formation of ices. 

The HDO/H2O and D2O/H2O observed ratios suggest that 
the ices present in IRAS 16293 and desorbing in the hot corino 
were made in a warm (~ 15 - 20 K) and dense (~ 10* - 5 x 
10* cm"^) translucent cloud. Recently, Taquet et al.[ ( |2013| mod- 
eled the formation of deuterated water ice in cold conditions 
with a pseudo time-dependent multilayer approach (see Taquet] 
|etal. [2013, for a full description of the model). They showed that 
the relatively high deuterium fractionation of water could be at- 
tributed to the formation of water ices in dark phases (either at 
hh ~ 10* cm"-* at low temperatures, ~ 10 K, and high visual 
extinction, larger than 4 magnitudes, or for higher densities, ~ 
10^ cm"-', at higher temperatures), where the D/H ratio starts to 
efficiently increase. Although the chemical models diverge on 
the deuteration processes which are enhanced for warm dust in 
one case (Cazaux et al.|201 l| l and for cold dust in the other case 
(Taquet et al. 20131, our results agree that the high deuterium 
fractionation originates from a dense medium (both models need 
dense medium for high D/H ratio). 

In the extended absorbing layer, which is possibly due to 
photodesorption mechanisms, the observed D2O/H2O ratio ap- 
pears much larger compared to the hot corino value, in contrast 
to the HDO/H2O ratio. It is also higher than all the chemical 
predictions. If the D2O/H2O ratio observed in the gas phase of 
this layer really reflects the D2O/H2O ratio at the grain surface, it 
would consequently suggest that water formed in a denser (> 5 x 
10* cm"-') region with warm (~ 20 K) gas. A high densit y (10^- 
10^ cm"^) is also required for the chemical modeling by Taquet| 
| et al.[ ( [20I3[ ). However, these results disagree with our obser- 
vational constraints obtained from the density of this absorbing 
layer (< 10^ cm"^, see Sect.pli. Nevertheless, we cannot exclude 
that, in this external layer, water can also be formed by gas phase 
mechanisms through ion-molecule reactions. This could also be 
one of the possible explanations of the difference between the 
D2O/H2O ratios determined in the hot corino and in the absorb- 
ing layer Other reasons can be mentioned to explain the dis- 
crepancy of the water deuterium fractionation between the hot 
corino and the absorbing layer During the heating of the grain 
surfaces, an isotope exchange between water molecules may also 
occur at 150 K ( [Smith et al.[[T997l [DuEeu et al.[[20T0| , leading 
to a decrease of the D2O abundance in the hot corino region. 
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Fig. 4. Deuterium fractionation of water predicted by the grain surface chemical model (see text) assuming 3 initial dust temperatures 
Jdust = 12 K (dashed lines), 15 K (dotted lines) and 20 K (solid lines) and densities n = 10^ (lower panel), 10"* (middle panel) and 5 
X 10^ cm"^ (upper panel). The HDO/H2O and D2O/H2O ratios determined at 3cr in the hot corino are indicated in pink with a red 
line for the best-fit value. The ratios of the absorbing layer are shown with a blue hne. 



Also, in the hot corino, large columns of H2O, HDO, and D2O 
could self-shield these molecules. Since H2O is more abundant 
than its deuterated forms, it should self-shield first, reducing its 
photodissociation. Then if important columns of HDO and D2O 
can be reached, these molecules will be shielded. This implies a 
stratification of molecules and should favor the amount of water 
compared to its deuterated forms in environments where lots of 
water is present in the gas phase. 



Recently, Persson et al. ( 2013| l used interferometric (SMA 
and ALMA) observations of one HDO transition (226 GHz) and 
two H2'^0 transitions (203 and 692 GHz) to determine the water 
deuterium fractionation in the hot corino of IRAS 16293. Water 
is only emitted by the core A and not by its binary component 
B. With their different method, jPersson et al. (2013) derive a 
fractionation ratio about a factor of 10 lower than ours ( |Coutens| 
|et al.|20l"2| . The origin of the divergence between these results 
is not clear, but could arise from the uncertain inner structure of 
the source. Indeed, a recent publication by Pineda et al. ( 20I2| l 
mentions that ALMA observations towards the source A would 
be consistent with the rotation of a disk, which is close to being 
edge-on. None of the two studies assume such a complex struc- 
ture in their modeling, maybe implying the divergence between 
the estimated HDO/H2O ratios. Compared with the predictions 



of the chemical model b y.Cazaux et al. (201 1 1, the ratio deter- 
mined by [Persson et al. ( 2013| l would require a low initial dust 
temperature (around 12 K), in contrast to our results. But doubly 
deuterated water has never been detected with interferometers, 
making it not possible to check whether the D2O/H2O ratio is 
also in agreement with an origin on cold dust grains and pro- 
vide more constraints on the initial density. Consequently, com- 
plementary D2O observations at high spatial resolution would 
be necessary to assess precisely the conditions in which ices 
are forming with a method based on interferometric data as in 
[Persson etaL] ( |20T3] l. 



5. Conclusions 



Thanks to the numerous HDO and H2'**0 transitions observed 
with the HIFI instrument on board the Herschel observatory, we 
were able ( [Coutens et al.[[2012) to constrain the abundances of 
these species from the hot corino throughout the envelope of 
the IRAS 16293 protostar to the external absorbing layer. In the 
present study, we have used the same non-LTE radiative transfer 
modeling to constrain the D2O abundances, employing the pre- 
vious detections of the para and ortho fundamental transitions 
( [Butner et aL]|2007[ [Vastel et aL[[2010| l as well as the tentative 
detection of the higher energy para-D20 2i,2-lo,i transition at 
898 GHz in Herschel/HIFI band 3b. To estimate the abundances, 
we used the full set of D2O collisional rate coefficients com- 
puted with ortho- and para-H2 ( Faure et al.|2012[ ), which were 
not available for the original detections. The results are presented 
in Table 2. The best-fit inner abundance Xin(D20) is about 7 x 
10"'", whereas the best-fit outer abundance XouiiJ^iO) is about 5 
X 10 '^. The absorbing D2O components cannot be reproduced 
without taking into account the external absorbing layer The 
column density in this layer is about 2.5 x lO'^ cm"^, leading 
to an abundance ranging from 6.6 x 10"'*' to 2.7 x 10"^ if the 
lines are produced in a photodesorption layer at Ay ~ 1^ magni- 
tudes. We then compared the HDO/H2O and D2O/H2O observed 
ratios with the predictions of a grain surface chemical model 
( [Cazaux et al.[[20lT| ), in which we considered the formation of 
ices by accreting species from the gas phase. Our estimations of 
the HDO/H2O and D2O/H2O ratios suggest that the ices present 
in IRAS 16293 and desorbing in the hot corino were made in a 
warm (~ 15 - 20 K) and dense (10"* - 5 x lO** cm"-') translucent 
cloud. The observed D2O/H2O ratio appears much larger in the 
absorbing layer than in the hot corino, and we explored tentative 
explanations. 
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